Wigner crystal versus Fermionization for one-dimensional Hubbard models with and 

without long-range interactions 
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The ground state properties of Hubbard model with or without long-range interactions in the 
regime with strongly repulsive on-site interaction are investigated by means of the exact diago- 
nalization method. We show that the appearance of A''-crests in the density profile of a trapped 
N-fermion system is a natural result of "fermionization" between antiparallel-spin fermions in the 
strongly repulsive limit and can not be taken as the only signature of Wigner crystal phase, as the 
static structure factor does not show any signature of crystallization. On the contrary, both the 
density distribution and static structure factor of Hubbard model with strong long-range interac- 
tions display clear signature of Wigner crystal. Our results indicate the important role of long-range 
interaction in the formation of Wigner crystal. 

PACS numbers: 03.75.Ss, 67.85.Lm, 71.10.Fd 



I. INTRODUCTION 

The Hubbard model has been widely used as a minimal 
model to describe strongly interacting electrons in a solid 
and played a particularly important role in understand- 
ing the physics related to the quantum magnetism and 
high-temperature superconductivity. Recent experiment 
progress in trapping the quantum gas in optical lattice 
Hl-Q has renewed the interest in this basic model d, Q . 
In comparison with real materials in solid, the fermionic 
quantum gas trapped in optical lattice is a more ideal re- 
alization of the Hubbard model [1] , which provides a lab- 
oratory to simulate and study this basic quantum many- 
body system from different aspects. Particularly current 
experiment techniques also prompt intensive study of 
quantum many-body systems in reduced one-dimensional 
(ID) systems, among which important benchmarks in- 
clude experimental realization of the Tonks- Girardeau 
(TG) gases @, 0], and the Mott transition in ID tight 
tubes [8|. Furthermore, the recent experimental realiza- 
tion of two-component Fermi gases in a quasi- ID geome- 
try with tunable interaction strengths by Feshbach reso- 
nance provides a unique possibility to experimentally 
study the ID Hubbard model even in the strong inter- 
action limit. The good tunability of optical lattices en- 
ables us to experimentally elucidate the subtleties of ID 
quantum many-body systems, for example, the physical 
properties in the strongly interacting limit. Due to recent 
experimental progress on chromium atoms [lol . flTj and 
heteronuclear polar molecule produce d by stimulated Ra- 
man adiabatic passage technique [l^ Il3| . dipolar atomic 
systems with strong long-range dipole-dipole interaction 
(DDI) have also attracted intensive studies [l^ - [l9| . In 
addition, schemes of creating Coulomb-like interaction 
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in cold atom systems have also been proposed [20|, l21| . 

In the presence of long-range repulsive interaction, an 
interesting issue for the ID fermionic system is the ex- 
istence of Wigner crystal phase For a ID gas of 
electrons with long-range Coulomb interaction which has 
been investigated by Schulz in the scheme of bosonization 
p^ . the extremely slow decay of density correlations at 
wave vector Akp has been used as a signature of appear- 
ance of Wigner crystal. Here kp ~ Trn is the Fermi wave 
vector with n the particle density. After that, explor- 
ing Wigner crystal phase in the short-range interacting 
fermion systems, for example the Hubbard model, has 
also been carried out |23l - l26l | . By using bosonization and 
density matrix renormalization group methods to study 
the one-dimensional finite Hubbard chains with hard- wall 
boundaries, a crossover from 2kF- to 4fci?-density oscil- 
lations with the increasing of repulsion has been found 
[lEIlSl: and sometimes the Akp density oscillations was 
taken as a signature of the emergency of Wigner crys- 
tal. In the limit of infinitely strong interaction, the on- 
site repulsion imposes an effective Pauli principle between 
fermions with different spins. Consequently, the ground 
state in the strong repulsion limit is degenerate for states 
with different spin configuration |27H29| . and thus the 
ground state density distribution of an equally mixed 
trapped system is identical to the fully polarized system 
with the appearance of the doubling of the crests in den- 
sity profiles. Such a phenomenon, which is also refereed 
as fermionization, was recently observed by the experi- 
ment in a two-particle system with tunable interaction 
using two fermionic ^Li atoms (soj . 

As the density distribution is not solely a convincing 
criteria of Wigner crystal, the density-density correla- 
tion function plays an im por tant role in characterizing 
the Wigner crystal phase [3l|. In the present work, we 
shall study both the density distribution and density cor- 
relation function of the Hubbard model either without 
or with long-range interactions. For the Hubbard model 
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with only the short-range interaction, we find that the 
density correlation function does not show any signature 
of Wigner-crystal correlations in spite of the fact that 
the ground-state density profile displays pronounced Akp 
Wigner oscillations in the strong repulsion limit. To un- 
veil the role of the interaction range, we then study the 
Hubbard model with additional long-range interactions 
and compare the calculated results with those of the Hub- 
bard model. We find that the system with long-range in- 
teractions indeed exhibits quite different behaviors from 
those of short-range Hubbard model. Particularly, the 
clear evidence in static structure factor has indicated the 
important role of the long-range interaction in the for- 
mation of Winger crystal. 

The remain of the paper is organized as follows. In 
Sec.H, we introduce our studied models and main con- 
clusions based on the calculation of models. In Sec. HI, 
we present our calculation results of density distributions 
and static structure factors in subsection A and B for 
the Hubbard model and Hubbard model with additional 
long-range interactions, respectively. A summary is given 
in the last section. 



II. MODEL HAMILTONIAN 

The basic model describing the short-range interacting 
fermions in a ID lattice is the Hubbard model [s^l 

H = ^(Cq+i^ + H.c.) + UJ2 ^»T"»i' (1) 

i(7 i 

where c|^(cj) is the creation(annihilation) operator of the 
fermion with the spin a, hi^ {nn) is the spin- up (down) 
particle number operator, t is the hopping strength which 
is set to the unit, and U is the on-site interaction 
strength. The ID Hubbard model can be realized by 
trapping a two-component Fermi gas in a deep ID opti- 
cal lattice. 

To study the long-range interaction, we consider the 
following Hamiltonian 

H = y^(cLc»+i<T + H.c.) + fii-fhii 

i(T i 

with V the strength of the long-range interaction and 
a the decaying exponent. For dipolar atoms confined 
in a ID optical lattices, the dipole-dipole interaction de- 
cays as l/|i — with the exponent a = 3 whereas 
a = 1 corresponds to the electron system with long-range 
Coulomb interaction. In the present work, we consider 
the case with repulsive on-site interaction U > and 
repulsive long-range interaction V > 0. 

In the following section, we shall first calculate the 
density distribution and the density-density correlation 



function for the Hubbard model ([T]) with the on-site in- 
teraction strength varying from weak to strong repul- 
sion limit by exact diagonalization method. In order to 
see the oscillation of the density profile, we shall use 
open boundary condition (OBC) in the calculation of 
density distribution, whereas periodic boundary condi- 
tion (PBC) shall be used in the calculation of structure 
factor. In the limit oi U — !• cxd, the total ground-state 
density distribution with different spin configurations is 
identical to that of the polarized noninteracting fermions, 
despite the fact that the spin-dependent densities are dif- 
ferent. This is consistent with the exact result obtained 
via the exact construction of many-body wavefunction 
[2^. Our numerical diagonalization results also show 
that the structure factor, which is the Fourier transfor- 
mation of density-density correlation function, in the in- 
finite repulsion limit also displays the same behavior as 
the fully polarized fermions. These results suggest that 
the appearance of the doubling of peaks of the equally 
mixing system (or Akp Wigner oscillation) in the infi- 
nite repulsion limit can not be taken as a hallmark of 
the emergency of Wigner crystal. It is just a kind of 
"fermionizaiton" between fermions with different spins 
as a result of the effective Pauli principle imposed by the 
infinite repulsion. 

The absence of Wigner crystal phase in the short- 
range Hubbard model implies that the long-range inter- 
action may play a very important role in the formation 
of Wigner crystal phase. To unveil its role, we then cal- 
culate the density distribution and the density-density 
correlation function for the Hubbard model with the ad- 
ditional long-range interaction described by Hamiltonian 
For conveniences, we shall focus our calculation to 
the case with a = 1 for which the long-range interac- 
tion decays more slowly than the dipole-dipole interac- 
tion. Keeping the on-site interaction fixed, we find that 
the Wigner crystal phase emerges with the increase of the 
long-range interaction strength V. Finally, we also calcu- 
late the case with a = 3 which exhibits similar behavior 
as the case of a = 1. 



III. RESULTS AND DISCUSSIONS 

A. Fermi-Hubbard model 

We first calculate the ground state density distribution 
ni = rii-f- -|- nil Fermi-Hubbard model ([1]) under the 

OBC, where Uia- — {fiia) is the spin-dependent ground 
state density distribution. In Figll] we show the density 
profiles n.i and rii^ versus different U for the balanced 
Fermi-Hubbard model ((J) with = N^^ = and the 
lattice size L = 24. For the balanced case, we always 
have = n^j^ = ?^i/2. In the limit oi U —^0, we have 

riia = Ym" \4>i?^ where 4>i{i) = ^ is the l- 
th single particle wave function of the open tight binding 
chain [33, which fulfills the OBC of <?ii(0) = 0i(i-fl) = 0. 
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FIG. 1: (Color online) The ground-state density distribution 
of the Fermi-Hubbard model for different U with A'^ = A'4, = 3 
confined in a lattice with L = 24 under the open boundary 
condition. Here, solid square is for Ui = riif + nii, solid circle 
is for and solid triangle is for riii . The density distribution 
jii in the limit of f/ 00 is shown in (d) marked by hollow 
inverted triangles. 



FIG. 2: (Color online) The density distribution of the Hub- 
bard model for different U with Nf — 2, A^'^ — 4 confined in 
a lattice with L = 24 under the open boundary condition. 
Here, square is for m, circle is for riif and triangle is for n^j,. 
The density rii for i7 — 00 is shown in (d) marked by hollow 
inverted triangles. 



It is obvious that there are three crests in the density dis- 
tribution of Hi for iVj- = Ni ~ 3 in the weakly interacting 
limit. As the interaction strength U keeps increasing, the 
density distribution displays clearly six crests for U = 10. 
In the limit of U ^ 00, we have ni = '^f' \(f>i{i)\'^ with 
N = Nf + Nf for different spin configurations accord- 
ing to exact construction of many-body wavefunctions 
[29| due to the induced Pauli exclusion principle between 
antiparallcl spins. As shown in FigHJ the total density 
profile for U ~ 200 obtained by the exact diagonaliztion 
method is completely overlapped with the exact result in 
the strongly interacting limit. For the imbalanced case 
with Nf = 2, Ni = 4 and L = 24, we show the change 
of density profiles of the Hubbard model with the in- 
crease of U in FiglH Despite the fact that now we have 
ni-f 7^ nil, total density distribution also displays 
clearly six crests for U = 10. Also, the total density 
profile for U = 200 coincides with the density profile of 
the fully polarized system. Our numerical results clearly 
show that the total density distributions of the interact- 
ing fermion systems with different spin co nfig urations are 
identical in the limit of infinite repulsion [29| . 

Next we calculate the static structure factor, which 
is defined as the Fourier transformation of the density- 
density correlation function: 

Sik) = ^5]e^'^-('-^)[(n,n,) - (3) 
where k = 2Tmr/L with m — 0,1,..., i. In general, the 



Wigner crystal phase can be characterized by the en- 
hancement of 4:kp peak in the the static structure factor. 
In Figl31[a), we show the change of the static structure 
factor with the increase of interaction strength U for the 
case of TV^ = TV; = 3 and L = 24 under the PEG. The 
imbalanced case of iV^ = 2 and A'^i = 4 is shown in 
FigEfb). For both cases, as U is large enough (for ex- 
ample, U = 200), the distribution of the static struc- 
ture factor approaches to the distribution function of a 
fully polarized TV-particle free fermion system as shown 
in FiglSJc). Except of the drops at kL /[2tt) = N,L-N, 
no obvious peaks in the the static structure factor are 
detected. Our numerical results clearly indicate that the 
appearance of N peaks in the ground state density pro- 
file of a finite trapped Hubbard system is not a mark of 
the emergence of Wigner crystal in the strongly repulsive 
limit. 



B. Hubbard model with long-range interaction 

Now we study the Hamiltonian of the Hubbard model 
with long-range interaction ([2]). For the case of y = 0, 
the ground state properties of the system almost does not 
change with the increase of on-site interaction strength 
when the U is large enough. In order to explore the ef- 
fect of the long-range interaction, we shall keep U fixed 
by taking U = 200 and change V in the following calcu- 
lation. For simplicity, we focus our discussion on the case 
of a = 1, and finally discuss the case of a = 3 briefly. 
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FIG. 4: (Color online) The density distribution of Hubbard 
model with long-range Coulomb interaction for different V 
with U = 200, A'^t = = 3, and L = 24 under the open 
boundary condition. 



FIG. 3: (Color online) Static structure factor S'(fc) vs kL/{2n) 
for different U with the periodic boundary condition, (a) is 
for N-f = = 2,,L = 24, (b) is the case of A^-^ =2,N^ = 
4, L = 24, and (c) is the case of polarized free fermion system 
with = 6 and L = 24. 



In FigHl we show the density distribution for the long- 
range Hubbard model with a = 1, N-^ = = 3 and 
L = 24 under the OBC. AtV = 0, six crests have already 
appeared for U = 200. Nevertheless, with the increase of 
V, the particles moves apart away from each other as far 
as possible to minimizing the total energy. As shown in 
the FigUl very sharp peaks have emerged for = 50 and 
the height of peaks increases with the increase of V. In 
contrast to the case oiV = with only small oscillations, 
the density profile of the system with a large V resembles 
six well-separated localized wave packets which distribute 
uniformly. 

To characterize the Wigner crystal phase, we fur- 
ther calculate the static structure factor of the system 



with L 



24, Nf 



3 under the PBC. As 



shown in FiglSj^a), obvious peaks have already emerged 
at fcL/(27r) = N, L - N ioi V = 5. As V keeps increas- 
ing, the height of peaks increases further and new peak 
emerges at kL /{2tt) = 2N . In FigOb), we calculate 
the static structure factor for a system with L = 32, 



3 under the PBC. Similarly, as the 



long-range interaction keeps increasing, a series of peaks 
emerge at fcL/(27r) = mN with m the integer. We also 
demonstrate data for the static structure factor S{k) ver- 
sus k/{mT) in FigE] for systems with L = 24, U = 200, 
V = 150 and TV = 2, 3, 4, 5, 6, where n = N/L is fill- 
ing factor of the system. As shown in FiglSl systems 
with different filling factors display similar behaviors in 
the regime with strong long-range repulsive interaction, 
i.e., peaks emerging at k/irnv) = 2m. For the different 
systems, the positions of the peaks only depend on the 
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FIG. 5: (Color online) Static structure factor S{k) vs kL/2-n 
for systems with U = 200 and different V under the periodic 
boundary condition, (a) is for L = 24, A^^ = A'j, = 3. (b) is 
for L = 32, = 1, A'; = 3. 



density n. 

Next we consider the Hubbard model with long-range 
dipole-dipole interactions described by the Hamiltonian 
of Eq[2] with a = 3. In FiglTl^a), we plot the density 
profiles of the Hubbard model with dipole-dipole inter- 



actions for a lattice system with L — 24, iV+ 



under OBC. The static structure factors of the same sys- 
tem with PBC arc plotted in FiglTJb). It is obvious that 
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FIG. 6: (Color online) Static structure factor S{k) vs k/{nn) 
(n = N/L) for systems with U = 200 ,V = 150, L = 24 and 
different A*' under the periodic boundary condition. 



both the density profiles and static structure factors dis- 
play very similar behaviors as these of the system with 
a = 1. A quantitative difference is that the height of 
peaks in FiglTl^a) and (b) is lower than that of the cor- 
responding system with Coulomb interaction shown in 
FigUand FiglSJa), as the long-range dipole-dipolc inter- 
action decays much faster than the Coulomb interaction. 
Despite the minor differences, our calculated results in- 
dicate the existence of the Wigner crystal phase for the 
Hubbard model with strong diploe-diploe interactions. 

As a comparison, we also study the extended Hubbard 
model described by 



cl^Ci+icr +H.c.)+C/ fii^nii+V fiihi+i, 



(4) 

where only the nearest neighbor short-range interaction 
is considered. Density profiles and static structure fac- 
tors of the extended Hubbard model with L = 24, 
= N^^ = 3 are shown in FigjTlJc) and (d), respectively. 
For F = 50, we observe that the oscillations in the den- 
sity distribution of the extended Hubbard model are en- 
hanced, but the wave packets are still overlapping, which 
is quite different from the density profile of the dipolar 
system with well-separated localized wave-packets. With 
further increasing the nearest neighbor interaction, no 
obvious change is found. The static structure factors 
also exhibit different behaviors from the long-range in- 
teracting system. As shown in Fig|7Jd), only small peaks 
emerge at fcL/2(7r) = A^, L — iV for = 50 and the height 
of peaks does not change obviously with further increas- 
ing V. The absence of peaks at fcL/(27r) = 2A^ implies 
that the crystal phase is not detectable. The different 
behaviors exhibited in systems with the short-range and 
long-range interactions indicate the important role of the 
long-range interaction in the formation of Wigner crystal 
phase. 
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FIG. 7: (Color online) (a) Density profile of the Hubbard 
model with dipole-dipole interaction under the open bound- 
ary condition; (b)static structure factor of the Hubbard model 
with dipole-dipole interaction under the periodic boundary 
condition, (c) density profile of the extended Hubbard model 
with open boundary condition; (d) static structure factor of 
the extended Hubbard model with periodic boundary condi- 
tion; Here, U = 200, L = 24, Nf = ^ 3. 



IV. SUMMARY 

In summary, we study the ground state properties of 
finite Hubbard systems either with or without long-range 
interactions by using the exact diagonalization method. 
We find that the appearance of A^-crests in the density 
profile of Hubbard model in the strongly repulsive limit 
can not be taken as the only signature of Wigner crys- 
tal phase, as it is induced by the effective Pauli principle 
between antiparallel-spin fermions enforced by the infi- 
nite repulsion. The absence of Wigner crystal phase is 
clearly verified by the calculation of the static structure 
factor, which shows the same behavior as the polarized 
free fermion system with no Akp peaks. If the long-range 
interaction is considered, the Hubbard model with long- 
range interaction can form a perfect crystalline phase in 
the regime with strong long-range repulsions. The exis- 
tence of the Wigner crystal phase can be well character- 
ized by its density distribution and the emergence of a 
series of peaks in the static structure factor. Our study 
unveils the important roles of the long-range interactions 
in the formation of the Wigner crystal phase. In view 
of the rapid progress in manipulating the dipolar atomic 
systems with tunable interactions, our results are possi- 
ble to be tested by the potential cold atom experiments. 
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